Superconductivity from Flat Dispersion Designed in Doped Mott Insulators 
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Routes to enhance superconducting instability are explored for doped Mott insulators. With the 
help of insights for criticalities of metal-insulator transitions, geometrical design of lattice structure 
is proposed to control the instability. A guideline is to explicitly make flat band dispersions near 
the Fermi level without suppressing two-particle channels. In a one- dimensional model, numerical 
studies show that our prescription with finite-ranged hoppings realizes large enhancement of spin- 
gap and pairing dominant regions. We also propose several multi-band systems, where the pairing 
is driven by intersite Coulomb repulsion. 
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Basic properties of doped Mott insulators have been a 
subject of recent continued studies One of the goals 
of the studies is to find ways to design instabilities such as 
magnetism and superconductivity by controlling material 
parameters in a realizable way. It is desired to control the 
instabilities by utilizing inherent character of the doped 
Mott insulators and the critical nature of metal-to-Mott 
insulator transitions. In this letter, possible prescriptions 
to control the instabilities are proposed. 

For single-band Hubbard and t-J models on a square 
lattice, simple scaling properties have been observed 
numerically near the transition from metals to Mott 
insulators implying nontrivial and singularly 

momentum-dependent correlation effects, where single- 
particle excitations around some particular points as 
(tt, 0) and (0, tt) in the momentum space play crucial roles 
with the emergence of a flat dispersion Q . This criticality 
extends to 20-30% doping range while appears only below 
some fraction of the scale of the exchange interaction J 
which is also the scale of the effective bandwidth. The flat 
dispersion around (7r,0) appears at similar scales. The 
flat dispersion was also numerically observed in the 2D 
Hubbard model in earlier works by Dagotto et al , Bu- 
lut et al. Q and Preuss et al. with comparison to uni- 
versally observed flat dispersions in angle-resolved pho- 
toemission data of the high-Tc cuprates . Although 
the observed flat dispersion is certainly resulted from a 
strong correlation effect with strong damping ||,|l^,^,0, 
the microscopic mechanism for the flat dispersion still 
waits for more complete understanding. Furthermore, 
its persistence in multi-band or more complex systems 
has not been well studied. 

In this letter we discuss that promotion of the above 
scaling behavior and the flat dispersion offers a way to 
control potential instabilities. We show that even when 
a flat hand dispersion is designed near the Fermi level 
by controlling lattice geometry and parameters, it en- 
larges the critical region under the suppression of single- 
particle coherence in the proximity of the Mott insulator 
mentioned above, thereby enhances the instability. One 
might argue that if a flattened dispersion is designed, it 
simply makes the correlation effects relatively larger only 
through the change in the ratio to the effective band- 



width with the enhanced density of states. This is, how- 
ever, not the whole story on the verge of itinerant and 
correlation-induced localized states. The metallic exci- 
tations are determined from the coherent one near the 
Fermi level while the two-particle processes including the 
superexchange interaction J is rather determined mainly 
from a local, incoherent origin in the real space when the 
electron correlation is strong. It opens a possibility of 
enhancing the two-particle instability by suppressing the 
dispersion only near the Fermi level simultaneously keep- 
ing the amplitude of the two-particle processes large. By 
the suppression of only the single-particle coherence, two- 
particle processes work selectively and effectively since 
competition processes such as pair-breakings and damp- 
ing by the single-particle channel are suppressed. We 
flnd no particular dependence on dimensionality for this 
mechanism. This is, however, not possible in single-band 
models with nearest-neighbor hopping because the two- 
particle processes is not independent of the band dis- 
persion near the Fermi level. For example, J and the 
so-called three-site terms in the t-J models are simply 
scaled by the transfer t as /U in the strong correla- 
tion expansion. We discuss below how the independent 
control can be made in more complex systems. 

The flattened dispersion makes degenerate excitations 
and may cause various instabilities. In this context, we 
note that ferromagnetic instabilities by the flat band has 
been extensively studied |p^|-p^. However, we consider 
here only the cases with the singlet ground state at half 
flUing, where the ferromagnetic instability is suppressed. 

In ID t-J models, the phase diagram in the param- 
eter space oi J /t and electron concentration n shows a 
general tendency of stronger pairing instability namely 
the larger Tomonaga-Luttinger exponent Kp with spin 
gapped excitation for larger J/t, though the phase sep- 
aration interrupts the enhancement |li^-[l5[|. Note that 
the pairing correlation is the most dominant if Kp > 1. 
The general tendency for enhanced pairing instability for 
large J/ 1 also holds in 2D . These are consistent with 
the mechanism we discussed above. We deflne the ordi- 
nary t-J model with the three-site terms: 
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where (ij) and (ijk) are nearest neighbors and cj^icia) 
creates(annihilates) electrons with constraint to exclude 
the double occupancy. If the band dispersion is sup- 
pressed only near the Fermi level, the effective transfer 
decreases while J can be retained mediated by incoher- 
ent but dispersive part. This is indeed possible by in- 
troducing longer-ranged transfer. We introduce the ID 
model with the third, fifth and seventh neighbor trans- 
fers, is, ts and t^, respectively and for the exchange part, 
Ji = JiU/t)"^ accordingly, in addition to t and J [gO|. We 
design flat nonintcracting dispersion e{q) by optimizing 
tsjts and to make terms up to the sixth order vanish 
in the wavenumber q around q = ±7r/2, the Fermi level 
at half filling. The spin-gap boundary and the exponent 
Kp = ^ttDu'^k/A were calculated by exact diagonaliza- 
tion of the Hamiltonian up to 16 sites, where the Drude 
weight D and the compressibility k were calculated at 
filling n following the procedure in the literature p3|-p^ . 
For the spin gap boundary, we used the level crossing 
method for accurate estimates [|l5|. The results have 
practically no system size dependence implying a reli- 
able estimate of the thermodynamic results. They show 
remarkable enhancement of both the pairing correlations 
and spin gap region as in Fig. |^, where the phase sepa- 
ration is absent. It shows the mechanism we proposed is 
indeed effective. We have also examined other types of 
band modifications and confirmed that the above mecha- 
nism mainly determines the enhancements pO| . We have 
found that the corresponding Hubbard-type model with 
the same distant-ranged transfers also show a similar en- 
hancement. However, the system sizes we could study is 
not large enough for reliable estimates of the thermody- 
namic limit. We note that the enhancement of pairing 
was also reported in the ladder model |r^-[l9[|, when the 
Fermi level lies near the top (bottom) of bands, where a 
flattening is present. Our flnding suggests one prescrip- 
tion, namely, look for a tuned flatter dispersion arranged 
only near the Fermi level of the doped Mott insulator in 
quasi one-dimensional conductors. Since the mechanism 
itself is not conflned to quasi-one dimensionality, it opens 
various possibilities with rich physics along this line. 

Next, we analyze several multi-band models with non- 
trivial band-flattening effects to get further insight on 
this issue and to provide hints for material synthesis. We 
introduce an extended Hubbard model 



Hh = -tY,{clcj, + H.c.) ~ t'Y^iclc^^ + H.c.) 

{ij) (Im) 

+ uJ2(^ii - hi^a -h + Yl ^v^i^j (2) 



where cl^{cia) represents the creation (annihilation) op- 
erator of an electron at site i with spin a and riio- = c|g.Cio- 
and rii = number operators. The lattices are 

defined below in Figs. § - 1| where t connects bonds with 
solid lines and t' connects broken bonds. 



The first example belongs to a category of one-quarter 
depleted square lattice. Our model is described by the 
Hamiltonian (2) with the lattice structure illustrated in 
Fig. ^. When t' = 0, it reduces to a lattice consid- 
ered by Lieb jl^ if isolated spins on a quarter of lattice 
points are depleted. It has ferromagnetic ground state 
at half filling as Lieb proved. For t' ^ 0, however, the 
ground state becomes singlet at half filling due to the 
Lieb-Mattis theorem ||2l[]. For smaller t'/i, the nonin- 
tcracting dispersion shows flattening in the middle two 
bands among four in total. The two flattened dispersions 



are given hy E = ±\ - + 2\/r)(S - 2Vr)), 



2u^(cos kr^ 
cos ky)^,u 



with ^ — (w^ + — u^) cos {kx 
cos ky) — — — w^, r = w^u^ 
2tt7VPTt^, u = Vt^ V = {f- t''^)/Vt^ + t'^. In 
the limit of t' 0, these two bands become completely 
flat. However flat the band becomes with small t' /t, the 
superexchange interaction on the solid bonds in Fig. | is 
unchanged and favors the antiferromagnetic order at half 
filling. This lattice structure may favor triplet pairing in- 
stability near t' = in contrast to the cases below, due 
to spin polarization of the fiat nonbonding band. 

The second model has slightly different lattice from the 
case above as illustrated in Fig. ^. In this bipartite struc- 
ture, the number of connected A and B sublattice points 
are equal even at t' — after removing isolated spins, 
where the ground state is singlet at half filling and the 
Mott insulating state with the antiferromagnetic order is 
expected in the thermodynamic limit. The flattening of 
the band is not complete in this case, but still has an 
extended region of flat plateau. 

The third example is illustrated in Fig. ^. The non- 
interacting bands consist of a dispersive bonding and a 
completely flat antibonding bands. Under electron dop- 
ing, carriers go into the flat band. Because of a strong 
frustration in contrast to the other cases, the magnetic 
correlation may be suppressed. 

For the fourth system, the lattice structure is illus- 
trated in Fig. p. The square lattice structure in Fig. ^ 
is not important and the following argument applies to 
any other bipartite lattices if the 3-site unit cell has the 
same structure. The band structure consists of anti- 
bonding, nonbonding and bonding bands from high to 
low energies given by Sa — —tCi + \J {tCiY + 2t''^, e„ — 
0, Sb = -tCi - yJ{tCiY + 2i'2, respectively with Ci = 
cos fcx-l-cos ky and an energy gap Ag = —21+^/4:1"^ + 21'"^. 
The nonintcracting ground state at half filling is given by 
filled bonding and half-filled nonbonding bands. 

Below we discuss a general aspect of interaction effects 
more or less valid in all the above cases although we first 
take the fourth system as an example. For nonzero C/, 
the flat band splits into upper and lower Hubbard bands. 
At half filling, the lower Hubbard band of this nonbond- 
ing band is filled leading to the Mott insulator. Each 
cell is occupied by precisely one electron in the nonbond- 
ing orbital and an exchange interaction between these 
nonbonding electrons in the higher order in U may sta- 
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bilize the antiferromagnetic order. In the perturbation 
expansion in [/, we have mixing of antibonding compo- 
nents into the filled bonding band. However, in the per- 
turbation, the number of nonbonding electrons at each 
cell can be changed only by two so that in all the cells, 
precisely one nonbonding electron is kept up to the in- 
finite order in U. Even away from half filling, when we 
once assign configuration for singly occupied and empty 
sites for nonbonding orbitals, they do not have dynamics, 
where macroscopic degeneracy remains. If the perturba- 
tion in U converges, the system remains insulating up to 
the doping concentration S = 1/3. Intracell but intersite 
interaction plays a similar role to U. 

An important process to kill the insulating state arises 
from the intercell Coulomb repulsion V|i_j| {i ^ j) be- 
tween the i-th and j-th cells. For simplicity, we take V 
between the same sublattice points 1, 2, or 3 in Fig. ||. 
Other combinations play a similar role. Away from half 
filling, there appear empty and doubly occupied cells in 
the nonbonding orbitals. In the first order perturbation 
in V|i_j| , a pair of electrons (holes) on the singly-occupied 
nonbonding orbitals at sites i and j each is excited to the 
itinerant antibonding (bonding) band, where the excita- 
tion energy is roughly 2Ag. In the second order in V , a 
pair hopping from i,j cells to Z,rn cells occurs if I and 
m cells have initially no nonbonding electrons. Here we 
show the perturbation expansion up to the second order 
for the third model (Fig. |): 
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s[k'{ri+s' - Tj+s)] cos[fc(n - n)] 
2t' - 4:t{cosk + cos k') 



where a\ ^ ^ creates a flatband electron at site i with 
spin a. Similar effective Hamiltonians are derived for 
the other models. 

From the second order process, the localized non- 
bonding electrons melt. It also induces superexchange 
type magnetic interaction from the charge diagonal part. 
When V has some extended range, the superconducting 
order with the pairing form factor with this range may 
occur. If the range becomes longer, the magnetic inter- 
action becomes more frustrated and the magnetic corre- 
lation is suppressed. The longer-ranged Coulomb inter- 
action also helps to inactivate the first order process in 
V . This is understood from the limit of infinite-range V . 
If does not depend on the combination the 

first order process just gives a constant independently of 
the electron configuration. If V\i-j\ is given by a screened 
form ~ [exp{—rij / R) / rij]Vo with the range i?, the sec- 
ond order process generates averaged kinetic energy in 
the order {Vi^ / J\g)5'^{R/ af[t/ i:^gf in 2D where a is the 
lattice constant. The average Coulomb repulsion energy 



generated from the first order process is roughly 
in d dimensions. Then the pairing interaction in the sec- 
ond order dominates for larger R. By a rough estimate 
inferred from the melting transition of the Wigner crys- 
tal j22|, the dominance may well be realized for R/a > 3 
a,t 6 ^ 0.1 and Vb/Ag ^ ^g/t ^ 1- If this condition 
is not well satisfied, the charge ordering may compete 
with the superconductivity. A related multi-band pair- 
ing mechanism was pointed out before pSlp^ ] , where the 
pair transfer through U drives the s-wave superconduc- 
tivity. Although several other proposals for pairing in 
multiband models with V have been made |^,^, they 
are rather different from ours in the sense that our mecha- 
nisms come from (interband) pair transfers from flattened 
band near the Fermi level to dispersive channel. An ar- 
tificial but helpful limit to understand the instability in 
our mechanism is the case of infinite-range transfer with 
the uniform amplitude t/^/N for system size N. In this 
case, the mean field solution with the order parameter 

A = j:kf(k)CkaC-ka' with /(fc) = E^Jy\^-3\e"'^'~'^/N 

becomes the correct description of the superconduct- 
ing ground state, where the Hamiltonian has the form 
H = —jjA^Ajj^. This gives an example of supercon- 
ducting order derived from purely repulsive interaction. 
Up to this discussion, the triplet and singlet pairings are 
degenerate. This is lifted by perturbation such as the 
frustrated magnetic interaction. 

It is tempting to say that the above designed flatten- 
ing seems to be an extreme limit of the high- Tc cuprates. 
In our models, the flat dispersion is separated from the 
dispersive bands. The high-Tc cuprates seem to be more 
similar to the fourth model at large U, where the disper- 
sive (anti)bonding and the flat nonbonding bands coexist 
at the Fermi level. 

It would be desired to seek for materials which follow 
the prescription presented here. For the flrst and second 
systems, regular 1/4 depletion of square lattice or substi- 
tution (for example one quarter of Cu ions with nonmag- 
netic ions on the Cu02 lattice) may satisfy the require- 
ment. It is interesting that this structure has some sim- 
ilarity to the charge ordered (or stripe) phase suggested 
in La-based high-Tc cuprates ||2^ . This connection pro- 
vides a new view on the stripe problem pTf . The fourth 
model may be realized by 3-layered structure with top 
and bottom layers designed to remove ligand atoms such 
as oxygen or substituted with different ligand atoms to 
suppress transfer through them in those planes. 

One might argue that, with flat dispersions, impurity 
or phonon effects were serious in real materials. How- 
ever, when the coherent single-particle excitations are 
suppressed, the localization effect is rather determined 
under the competition to the two-particle process. The 
localization effect is reduced if the two-particle hopping 
process is retained larger. Anyhow impurity and phonon 
effects near the Mott insulator in the present situation 
are important open problems left for future studies. 

Our proposal to enhance the pairing instability is sum- 
marized: Design flat band near the Mott insulator by re- 
taining the preferred two-particle process. Due to the en- 
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hanced degeneracy and the suppression of single-particle 
process, the doped system gets stronger instability in 
the two-particle process. We have numerically shown 
that this mechanism and prescription work in an exam- 
ple of ID models with longer-ranged transfers. We have 
proposed several multi-band lattices, where the intersite 
Coulomb repulsion generates such two-particle processes 
through interband pair hoppings. Because the coherent 
single-particle process is absent, small intersite interac- 
tion immediately corresponds to the strong coupling limit 
of the pairing where the Fermi liquid state does not exist. 
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FIG. 1. Phase boundary of the spin gap region (circles) and 
the contour lines (with squares) for the TomonagarLuttinger 
exponent Kp = 1. The case for the t-J model with the 3-site 
term (open symbols) and that of the same model but with 

optimized dispersions by t-^ = 0.6, ts — 0.2 and tr = 1/35 and 
resultant J3,J5 and Jr, where Ji = Jtl (filled symbols). 
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FIG. 3. The sec- 

2. (a)The first q^^^^ model with alternating 
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FIG. 4. The third model, ladder structure with diagonal 
transfer. 
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FIG. 5. The fourth model, decorated square lattice. 
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